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Abstract: Interest in the protein folding problem has motivated a wide range of theoretical 
and experimental studies of the kinetics of the collapse of flexible homopolymers. In this Paper a 
phenomenological model is proposed for the kinetics of the early stages of homopolymer collapse 
following a quench from temperatures above to below the 9 temperature. In the first stage, nascent 
droplets of the dense phase are formed, with little effect on the configurations of the bridges that join 
them. The droplets then grow by accreting monomers from the bridges, thus causing the bridges to 
stretch. During these two stages the overall dimensions of the chain decrease only weakly. Further 
growth of the droplets is accomplished by the shortening of the bridges, which causes the shrinking 
of the overall dimensions of the chain. The characteristic times of the three stages respectively scale 
as N°, iV 1 / 5 and N 6 ^ 5 , where N is the degree of polymerization of the chain. 

PACS: 61.25.Hq, 83.10.Nn 
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I. INTRODUCTION 

The purpose of this Paper is to discuss the evolution 
of the structure of a single polymer molecule after the 
environment has been changed from a good solvent to 
a poor solvent. Our approach is via a phenomenolog- 
ical model for the kinetics of the early stages of col- 
lapse. We have in mind a long flexible polymer chain 
immersed in a simple solvent, and a change in the tem- 
perature T as the means of altering the solvent quality. 
The initial state is one of equilibrium in a good solvent, 
in which the chain is swollen and has radius of gyration 
R = Rp ~ TV 3 / 5 , where TV is the degree of polymeriza- 
tion. The final state is one of equilibrium in a poor sol- 
vent, in which the chain is collapsed into a dense, spher- 
ical globule, i.e., R — R c ~ TV 1 / 3 . The variation of the 
equilibrium state with the solvent quality is, with some 
important caveats, well understood Jl|-^|]. In marked con- 
trast, there is no comparable consensus concerning the 
kinetics of collapse following an abrupt T quench, i.e., a 
rapid change in temperature from above to below the 9 
temperature. The experimental picture is not clear be- 
cause of difficulties due to the competition between the 
collapse of the individual chains and the aggregation of 
different ones. As a result, there have been only few ex- 
perimental studies of the kinetics of collapse of isolated 
chains . Theoretical studies have utilized a variety of 
approaches in focusing on different aspects of the prob- 
lem. Langevin models p|-^0|, phenomenological mod- 
els [|lT}-[l3f and computer simulations |]l4]-f20| have been 
used to consider the kinetics of collapse in the absence 
of topological constraints. Other studies have focused on 
the role of internal entanglements [ pTU2^| o r the compe- 
tition with chain-chain aggregation pap4| . In spite of 
the extensive activity is this area, a complete picture of 
the kinetics of collapse has yet to emerge. The search for 
such picture is motivated primarily by the intense current 
interest in the protein folding question p5| p7| because 
early stages of protein folding are thought to proceed as 
the collapse of flexible homopolymers. 

The collapse of a flexible chain following an abrupt 
change of solvent quality displays certain features of a 
first-order phase transition. These features occur for 
chains of finite TV for which the surface free energy of 
the globule plays an important role. Two such features 
are important to the present discussion. First, the onset 
of collapse involves, as we shall discuss, the formation 
of droplets of the dense "phase". Second, the collapse 
proceeds via intermediate states corresponding to the 
"tadpole" configuration of a stretched, collapsed glob- 
ule |2g], which involve the coexistence of dense globules 
and strongly stretched segments. With this picture in 
mind, it is possible to distinguish between four stages 
in the kinetics of collapse (Fig. [I]). Each stage is char- 
acterized by a length-scale and a time-scale. The char- 
acteristic time scale increases with the relevant length 
scale. The formation of p (<C TV) droplets of the dense 



phase, which we refer to as "pearls", is the fastest pro- 
cess. This "pearling" stage involves local rearrangement 
of the chain configurations within the randomly placed 
droplets during which the dimensions and configurations 
of the chain as a whole are only weakly modified. The 




FIG. 1. Proposed sequence of states during the collapse of 
a flexible homopolymer: 

(a) self-avoiding configuration under good solvent conditions; 

(b) pearling; 

(c) bridge stretching; 

(d) collapse of the pearl necklace; 

(e) equilibrium collapsed configuration in a poor solvent. 

nascent droplets then grow by accreting monomers from 
the bridges that connect them. The positions of the grow- 
ing droplets are roughly stationary and their growth is 
accommodated by the gradual stretching of the bridges. 
This "bridge-stretching" stage ends when the tension in 
the bridges attains the equilibrium value f co correspond- 
ing to the "tadpole" configuration of a stretched glob- 
ule. During this process the dimensions of the chain as 
a whole remain essentially constant i.e, R ~ TV 3 / 5 , but 
the configurations of the bridges are qualitatively mod- 
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ified. The number of droplets p is set by the require- 
ment that R ~ TV 3 / 5 at the end of this stage, thus lead- 
ing to p ~ iV 4 / 5 . Beyond this point, the growth of the 
"pearls" is accompanied by shortening of the bridges and 
the overall shrinking of the chain. The third stage, the 
collapse of the "pearl necklace" , involves the chain as a 
whole. The corresponding length scale iV 3 / 5 ) is 

larger and the characteristic time is thus longer. It is 
assumed that p remains roughly constant during these 
three stages. Eventually, the droplets come into contact, 
and coalesce into a single globule. Our discussion is lim- 
ited to the early stages of the collapse. The last regime 
will not be considered in the present paper. Within this 
picture, the characteristic times for the first three stages, 
pearling, bridge-stretching and the collapse of the pearl 
necklace are, respectively, r p ~ N°, Tb a ~ iV 1 / 5 and 
r m ~ N 6/5 - The model we propose differs from the 
existing phenomenological descriptions in two respects. 
First, the two initial stages were not considered in these 
models. Second, the treatment of the third stage, the col- 
lapse of the pearl necklace, is different. A more detailed 
comparison with the phenomenological models, as well 
as the Langevin-type theories and to simulation studies, 
will be presented in the final Section. 

The Paper is organized as follows. To obtain the tem- 
poral evolution of the various stages we balance the driv- 
ing force with the entropy production. A discussion of 
this approach is presented in Section II. In Section III 
we recall some key results concerning collapse blobs and 
the thermodynamics of collapse. In addition, we discuss 
the initial role of the collapse blobs as spinodal "modes" 
and the equilibrium deformation behavior of a collapsed 
chain. The collapse of a chain from a good solvent to a 
9 solvent is discussed in Section IV. The results are then 
used to analyze the first stage of the collapse, i.e., the for- 
mation of nascent droplets or "pearling" . In Section V 
we discuss the collapse of a chain with constrained ends. 
This discussion is then utilized to describe the second 
stage of collapse, in which the stationary droplets, pearls, 
grow by accreting monomers from the bridges thus caus- 
ing the bridges to stretch. The shrinking of the resulting 
string of pearls is analyzed in Section VI. 



II. SOME GENERAL CONSIDERATIONS; 
ENTROPY PRODUCTION METHOD 

The design of the phenomenological model presented 
in this paper, involves two steps. In the first step we pro- 
pose a sequence of stages followed by the collapsing chain 
(The stages involved were outlined in the introduction.) 
The second step is to calculate the characteristic time 
associated with the various stages in the proposed route. 
A simple method for calculating the characteristic times 
was proposed by de Gennes (yj . It focuses on a volume 
element of solvent containing a single chain and coupled 
to a thermal reservoir that maintains a temperature T. 



The systems as a whole, is closed. During collapse, the 
free energy of the chain F c hain decreases. The decrease 
in Fchain gives rise to an increase in Stotai the total en- 
tropy of the system as a whole. Since the collapse is the 
only irreversible process involved the changes in F c hain 
and Stotai are related by 

AF c hain = —TAStotal- (2.1) 

The dynamics of the process are described by 

dF c hain ^dStotal ,~ ~\ 

~^r~ = ~ ~dT- [ ' 

To implement this argument, it is necessary to specify: (i) 
Fchain as a function of the chain radius R. One assumes 
that dFchain/dt = (dF chain / dR)V, where V = dR/dt. 
The form of F c ham for the various stages will be dis- 
cussed at the appropriate sections. One must also specify 
(ii) the operative dissipation mechanisms that determine 
dStotai/dt in terms of R and V . This entropy production 
approach is independent of the model, in that its appli- 
cability is not limited to a specific choice of F c ham or of 
the dissipative modes. Nor does it provide guidance for 
choosing this crucial input. 

Two dissipation mechanisms are considered for the var- 
ious stages of collapse. First is the viscous dissipation 
associated with the motion of the solvent. Second is the 
frictional dissipation due to the Stokes drag on the chain. 
The magnitude of the viscous force (2{| acting on a unit 
volume is on the order of 77V 2 V where 77 is the shear vis- 
cosity. In the following we estimate it by rjV/R 2 . The as- 
sociated dissipation per unit time is thus roughly r/V 2 / R 2 
and the total hydrodynamic dissipation within a volume 
R 3 is TdSh as given by 

The Stokes drag force acting on a sphere of radius r mov- 
ing with a constant velocity V is 6irr]Vr, and the associ- 
ated "frictional" dissipation per particle is 

T^^vV 2 r. (2.4) 

To obtain the total frictional dissipation it is necessary to 
sum the contributions due to the hydrodynamically im- 
penetrable objects involved in the stage considered. The 
number of such objects, their radius r and the relation- 
ship between r and R differ in the various stages of col- 
lapse. The de Gennes argument is equivalent to equating 
the frictional force to the contraction force. It allows us 
however, to consider on equal footing the dissipation due 
to Stokes drag forces and to the hydrodynamic flows. Fi- 
nally, note that the dissipation increases with the length 
scale. This observation motivates the proposed sequence 
of stages which involve processes taking place at progres- 
sively larger length scales. 

The familiar methods of irreversible thermodynam- 
ics [|30) focus on linear responses, the Onsager relations, 
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and their refinements. In marked difference, the de 
Gennes approach is not limited to a linear regime nor 
does it relay on the analysis of the regression of fluc- 
tuations. The method applies to the complete dynamic 
range, provided that the driving forces and the dissipa- 
tion mechanisms are correctly identified. The dissipative 
modes invoked in our discussion hold under the assump- 
tion that the hydrodynamic description is applicable to 
the collapse of a polymer. It is expedient to use this 
approach, but one should note that the validity of con- 
tinuum theories at this length scale is a delicate issue. 



III. COLLAPSE BLOBS, SPINODAL 
DECOMPOSITION AND THE EXTENSION 
ELASTICITY OF GLOBULES 

Whilst the configurations of polymers in good and poor 
solvents are well understood, the nature of the transition 
between the two regimes is not fully resolved. Different 
models suggest that the transition may be continuous or 
discontinuous, depending on N and the stiffness of the 
chain. The issue of the order of the transition is, however, 
irrelevant when considering a collapse following a quench 
to the poor solvent regime. In this regime, the equilib- 
rium configuration of the chain is determined primarily 
by the Flory interaction free energy Fi nt . The elasticity 
of the chain, i.e., its configurational entropy, plays a rel- 
atively minor role, limited to determining the structure 
of the globule-solvent interface || . In the case of a single 
chain, for which the translational entropy of the polymer 
is irrelevant, the N — > oo limit of Fi nt determines the 
equilibrium behavior: 

^ = (l-0)ln(l-0) + X </»(l-0). (3.1) 

Here, 4> is the monomer volume-fraction and the Flory 
parameter, is a measure of the strength of the interac- 
tions giving rise to the mixing enthalpy. For a given 
polymer-solvent system, x depends only on the temper- 
ature T. In the good solvent regime, for which \ < 1/2, 
Fi n t is a convex function of <f>. A critical point occurs 
at (xc, 4>c) specified by Xc = 1/2 and <j) c — 0. When 
X > 1/2 the solvent is poor and the plot is concave in 
the range < <f> < <f>+. A common tangent construction 
between the origin and 0+ indicates the coexistence of a 
dense phase, of concentration cj)+ and the neat solvent. 
The 4>+ phase corresponds to the collapsed globule. This 
amounts to setting the osmotic pressure of the polymer 
7r within the dense phase to zero, as is expected because 
there are no free polymers in the solvent. The vanishing 
of 7r leads to ln(l — <f>) = <f> + x4> 2 - I n t ne vicinity of the 9 
temperature, the ln(l — <fi) term may be approximated by 
the first three terms of its Taylor expansion. The n = 
requirement then leads to 

<t>+ « -v/2w, (3.2) 



where v = a 3 (l — 2%) < is the second virial coefficient, 
a is the size of the monomer, and w = 1/6 > is the 
third virial coefficient. Thus, is determined by the 
balance of attractive binary monomer- monomer interac- 
tions and higher-order repulsive interactions. The chain 
elasticity does not play a role. In turn, <p + sets the ra- 
dius of the collapsed globule via Na 3 /R 3 sa <f>+, leading 
to R/a « (N/^+f/ 3 . The occurrence of an unstable re- 
gion in Fi nt (4>) is a crucial feature of this discussion. Due 
to this feature, the collapse of a quenched chain exhibits 
certain characteristics of first-order phase transition. 

The structure of the interface, and the corresponding 
surface tension, can not be obtained from the discussion 
presented above. A "blobological" argument can provide 
the missing information || . The configurations of a col- 
lapsed chain are characterized by a correlation length, £ c . 
On length scales smaller than £ c the monomer-monomer 
attraction is too weak to perturb the configurations of 
the chain, and it behaves as a random walk. This allows 
one to define collapse blobs, which comprise g c monomers 

1 /2 

such that g c a « £ c . A fully collapsed chain can be envi- 
sioned as a spherical globule composed of close-packed £ c - 
blobs. The correlation length £ c is determined by equat- 
ing the interaction free-energy due to monomer- monomer 
attraction in such a chain segment to kT. The interac- 
tion free-energy is vg 2 /^ 3 , where v is the second virial 
coefficient. In the vicinity of the 8 temperature, v may 
be expressed as v m a 3 (AT/ 6), where AT = (T - 6)/6, 
thus leading to 

9^{-^)\ W-(^). (3.3) 

The density ga 3 /£ 3 within a single blob is comparable to 
that of the dense phase <p + resulting from the phase sep- 
aration of free polymers in a poor solvent. The £ c -blobs 
attract each other, and the unperturbed chain forms a 
collapsed globule of radius R « (N / gc) 1 / 3 ^. A collapse 
blob in contact with the solvent is assigned an energy of 
kT. Thus, the boundary of the globule is associated with 
a surface free energy 

7 « kT/e c - (3.4) 

Note that the blob picture, as described above, is valid 
only when the chain comprises many collapse blobs, i.e., 
when the depth of the quench is such that 9 AT 
Q/N 1 / 2 and 1 < g c < N. We should also add that in 
the following we treat the £ c -blobs as hydrodynamically 
impenetrable spheres. 

When N is sufficiently large, the collapsed globule may 
be viewed as a droplet of the dense phase with density <fi+ 
in coexistence with the neat solvent. Within our model, 
the first stage of collapse involves 0+ droplets, pearls, 
that are joined by bridges which retain the configura- 
tions found in good solvent conditions. This situation is 
reminiscent of the condensation of a fluid droplets from 
a supersaturated vapor Q. In this last situation it is 
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found that droplets whose radius is smaller than a criti- 
cal radius r c tend to evaporate whereas larger drops, i.e., 
those for which r > r c , grow indefinitely. This similar- 
ity wrongly suggests that the formation of nascent pearls 
is analogous to the nucleation of critical droplets. This 
analogy is discredited by the following rough argument. 
The thermodynamic potential of a dense droplet coexist- 
ing with "free" monomers is given by 



A = F(N d ) - N df i, 



(3.5) 



where N d is the number of monomers in the droplet, 
and p, is the chemical potential of the free monomers. 
For a sufficiently large droplet F = NdUd + jr 2 , where 
fid is the chemical potential of a monomer within the 
droplet. Choosing the reference state to be a blob within 
the droplet and assigning kT to an "uncondensed" blob, 
leads to 



A_ 

kT 



Nd 
9c 



Nd 
9c 



2/3 



(3.6) 



which has a maximum at Nd ~ g c . As Eq. (3.6) is only 
valid for N d ^S> g c , this establishes only that there is 
no critical droplet of larger size. This argument does 
however suggest that chain segments incorporating g c 
monomers play a special role in the onset of collapse. 
This view is consistent with the definition of a £ c -blob: 
The monomer-monomer attraction is too weak to per- 
turb the chain configurations on smaller length scales, 
but is sufficient to induce collapse on length scales larger 
than £ c . It is also helpful to note that the shrinkage of a 
chain segment into a collapse blob does not encounter a 
free-energy barrier. Accordingly, a collapse blob can not 
be viewed as a critical droplet. Rather, the formation 
of a £ c -blob resembles spinodal decomposition with the 
caveat that the associated dissipation arises from con- 
figurational changes of the chain as a whole. Thus, in 
order to minimize the entropy production the first step 
of the collapse involves the formation of a small number 
p (<C N/g c ) of collapse blobs instead of a string of N/g c 
£ c -blobs 

The concave region in Fi nt is also the origin of the 
force-law governing the extension a collapsed globule. 
When the end-to-end distance of the stretched globule 
is fixed, a dense globule coexists with a stretched string 
of £ c -blobs. The tension f co associated with this co- 
existence is the driving force for the third stage of the 
collapse. The extension of the globule actually involves 
three regimes |2^j3^]. In the first regime the spherical 
form of the globule is initially deformed into a prolate 
ellipsoid, but maintains a constant volume correspond- 
ing to close packing of the £ c -blobs. Within this linear- 
response regime the free-energy penalty F incurred is due 
to the increase of the surface free energy 



F « jAA « 7 (L - R c ) 2 



(3.7) 



where A^4 w (L — R c ) 2 is the surface area increment 
associated with the deformation, and L is the length of 



the axis of rotation. The corresponding restoring force, 
/ = —dF/dR, is proportional to the strain (L — R c ) : 



f^- 7 (L-R c ). 



(3.8) 



This type of process can not proceed indefinitely. If it 
were pursued indefinitely, the distorted globule will as- 
sume a cylindrical shape and, eventually, form a string 
of £ c -blobs. This scenario gives rise to a van der Waals 
loop in the / vs. R diagram. In turn, this is indicative of 
instability with respect to a coexistence of a weakly elon- 
gated globule and a stretched string of £ c -blobs. This 
effect is reminiscent of the Rayleigh-Plateau instabil- 
ity f33]| involving the break-up of a fluid jet into a succes- 
sion of droplets. In the second, coexistence, regime the 
chain comprises a stretched string of n/g c £ c -blobs and 
a roughly spherical globule of (N — n)j g c closely packed 
£ c -blobs. To characterize this equilibrium situation it 
is necessary to minimized the free energy of a stretched 
chain comprised of a dense globule and a string of n/g c 
collapse blobs, i.e., 



F 
kf 



N -n 

gc 



N-n 
9c 



2/3 



T 2 

( 3 - 9 ) 



The first term allows for the transfer of the collapse blobs 
from the solvent into the "dense phase,' i.e., the interior 
of the globule. This is the primary driving term for the 
growth of the globule. The second term reflects the sur- 
face free energy of the globule. The elastic free energy of 
the string of blobs is described by the third term. Note 
that na 2 ~ (n/g c )^ 2 is the unperturbed span of an ideal 
string of collapse blobs. Making F stationary with re- 
spect to n, dF/dn = 0, leads to 



1 

9c 



N 



-1/3 



1 



L 2 



g c n 2 a 2 



(3.10) 



In the limit of large L and for n < N this leads to 
L n 

ic^ g c ' 



(3.11) 



The equilibrium free energy of the resulting tadpole con- 
figuration is, up to numerical prefactors 



F 
kf 



N-n 
9c 



2/3 



n N\ n 

+-. (3.12) 

9c g c J gc 



The first term allows for the surface free-energy yrz of a 
globule of radius r globule w [(AT - n)/5 c ] 1/3 £ c . The trans- 
fer free-energy of blobs into the dense phase is reflected 
in the second term. The third term is the stretching en- 
ergy of the extended string of length L s w (nj 'g c )£, c a s 
given by the Gaussian expression for the elastic penalty 
kTL 2 /na 2 . It is also possible to interpret this last term 
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as the surface energy of the extended string. The end- 
to-end distance reflects two contributions 



N-n 

9c 



1/3 



n 
9c 



(3.13) 



The first term is the radius of the globule, and the second 
is the span of the stretched string of £ c -blobs. When n 
is sufficiently large we may approximate dR as £, c dn/g c 
and the corresponding force law / = —dF/dL reads 



r 



globule ' 



(3.14) 



It is important to note that / decreases as r g i b u i e ap- 
proaches £ c . This finite-size correction is reminiscent of 
the Laplace law for the vapor pressure of small droplets. 
It is negligible over a wide range of extensions if N is 
large enough. Because of this correction different defor- 
mation scenarios are expected for the / = const' and the 
L = const' ensembles. In the first case the globule unrav- 
els completely as soon as a critical force f co /kT « l/£ c 
is applied. No globule-coil coexistence is expected. Such 
a coexistence is however expected when an end-to-end 
distance is imposed. In this case the tension in the 
string is f co . The onset of the coexistence occurs when 
L m R c + £ c as can be s een by equating f co /kT to the 
force given by Eq. ( |3.8D . The upper boundary of this 
second regime corresponds to a fully extended string of 
blobs, L max » (N / g c )^ c . Stronger extension is charac- 
terized by simple Gaussian elasticity, up to the onset of 
finite extensibility effects. 



to the dissipative losses associated with the shrinking, 
—TdStotai/dt, where dStotai/dt is the corresponding en- 
tropy production. Two dissipation processes are in- 
volved. One is the viscous dissipation due to the hy- 
drodynamic flow resulting from the contraction of the 
coil. As was discussed previously, the resulting entropy 
production within the coil is roughly 



(4.3) 



where V — dR/dt is the velocity of the flow. Second 
is the entropy production due to the Stokes drag force 
experienced by the monomers as they move through the 
solvent. The associated entropy production is approxi- 
mately 



T 



dS s 
dt 



nNaV 2 



(4.4) 




IV. COLLAPSE TO 9 CONDITIONS AND THE 
FORMATION OF NASCENT DROPLETS 



A quench from T > 6 to T 

,1/5 A r3/5 



causes the chain to 



shrink from R F / a w yVSjy 3 / 5 to R Q /a w N 1 / 2 (Fig. |). 
The dynamics of this process are of interest for two rea- 
sons. First, as we shall discuss, the formation of nascent 
pearls follows similar kinetics. Second, the discussion of 
this process is especially simple, thus providing a useful 
introduction to the subsequent sections. The quench to 
T = 9 turns off the repulsive interactions between the 
monomers so that driving force for the "collapse to 6>" 
is due solely to the Gaussian elastic free energy of the 
F e i or 



kT 



a 



In a, 



(4.1) 



where a = R/R and R = R(t) is the instantaneous 
chain radius at time t. To obtain the characteristic time 
for this process we equate the time derivative of the elas- 
tic energy, 



FIG. 2. Collapse to 6 conditions: 

(a) initial, self- avoiding configuration; 

(b) the final, random walk configuration. 

The dynamics of the collapse are determined by 
dF e i/dt = -T{dS h /dt + dS s /dt). Whilst in this case it 
is possible to solve the full equation, it is easier, and 
more instructive, to consider dF e i/dt = —TdSh/dt and 
dF e i/dt = -TdSs/dt separately. dF el /dt = -TdS h /dt 
leads to 



dt 

TZ 



or da 



(4.5) 



and the characteristic time associated with the hydrody- 
namic dissipation is thus the Zimm time of the ideal 
coil, i.e., 



yR 3 o 

kT 



(4.6) 



The frictional dissipation is dominant as it scales with N 
rather than R. It leads to 



dF el 
dt 



da 
~dt' 



(4.2) 



dt 



ada 



TZ 



(4.7) 
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and the corresponding characteristic time is the Rouse 
relaxation time JlJ of an ideal coil tr ps r]N 2 a 3 /kT or 



tr « t z N 



1/2 



iV 2 . 



(4.8) 



As the Rouse time is much longer, it dominates the pro- 
cess. 

Within our picture, the first stage of the collapse is the 
formation of nascent pearls comprising each of roughly g c 
monomers. The process involves p randomly distributed 
chain segments. Initially, each of the chain segments 
exhibits self-avoidance. Once the solvent is quenched, 
the monomer-monomer repulsions disappear and the seg- 
ments shrink to their ideal coil configuration. As the 
nascent pearls consist of only g c monomers the monomer- 
monomer attraction is not strong enough to perturb the 
ideal coil configuration within them. So long as p is small 
enough, the overall radius of the chain is only weakly 
modified and the dissipation is due to the local config- 
urational changes. In this picture it is possible to ne- 
glect the couplings between the droplets and between the 
droplets and the chain as whole. Thus, the dynamics of 
this process are essentially identical to those of the col- 
lapse of a chain comprising g c monomers to 9 conditions. 
The characteristic time is, accordingly the corresponding 
Rouse time or 



1/2 



2 

kT 9c 



AT 



(4.9) 



ends constrained to fixed positions separated by a dis- 
tance L ^S> R c (Fig. ||). The free energy of such a chain, 
neglecting the surface free-energy, is given by 



kT 



N -n 
9c 



l? 



(5.1) 



where n is the number of "free" monomers that have not 
been assimilated into a globule. We assume that the n 
monomers behave as a random walk, i.e., the instanta- 
neous span of the chain perpendicular to the end-to-end 
vector is 



R± 

and as a result 



(n/ 5c ) 1/2 £c 



n l ' 2 a, 



dn 



1 dR ± 



(5.2) 



(5.3) 



tension kT/t, 





where 



(a) 



H *~ 



(b) 



rtff (_8_\ 3 
kT \AT J ' 



(4.10) 



Note that r p is in- 



is the Zimm time of a collapse blob, 
dependent of N and p as expected for a local process. A 
short relaxation time scaling as N° is also predicted by 
Langevin type theories of collapse ||Eo| . 



V. COLLAPSE OF A CHAIN WITH 
CONSTRAINED ENDS: THE 
BRIDGE-STRETCHING REGIME 

In the second stage of the collapse, the roughly station- 
ary droplets, "pearls" , grow by "sucking in" monomers 
from the bridges joining them. This stage ends when 
the initially undeformed bridges assume the configura- 
tion of a stretched string of £ c -blobs. In comparison to 
the "spinodal" stage described in Section IV, this pro- 
cess is distinguished by two features. First, the driving 
force is now due to the transfer free energy of the free 
monomers into the dense phase. Second, the evolution of 
the droplets at this stage strongly affects the configura- 
tion of the bridges. The bridges stretch as the number of 
monomers in them decreases. We first consider the some- 
what simpler situation of the collapse of a chain with its 



FIG. 3. Collapse of chain with constrained ends: 

(a) the initial, swollen configuration; 

(b) the final "tadpole" configuration in a poor solvent. 



Relations (5.2) and (5.3) are actually invalid at the initial 
stages of this regime, when the configurations of the chain 
still exhibit self-avoidance. They are, however, expected 
to apply at the final stages, which dominate the relax- 
ation time. Invoking them implies that the relaxation of 
the chain configuration when n decreases is effectively in- 
stantaneous. The driving force for this stage, as specified 
by Eqs. (pD-(pl), is 



1 dF, 



chain 



kT dt 



L 2 



1 D dR x 
a 2 ± ~dT' 



(5.4) 



The driving force is opposed by the two dissipation mech- 
anism invoked in Section IV. In the present situation, 
however, the dissipation reflects the cylindrical symme- 
try of the system. The hydrodynamic dissipation, due to 
the flow of the solvent as it is expelled from the region 
occupied by the coil, is T(dS h /dt) ps rj(dR±/dt) 2 L. The 
Stokes drag on the bridges formed by the collapsed blobs 
is T{dS b /dt) ps Ti(n/g c )£ c (dR ± /dt) 2 . The dynamics of the 

collapse are determined by dF c h a in/dt = — T(dS 'total /dt). 

l li 

At equilibrium, n is n eq ps (L/^ c )g c pa g c (L/a). Intro- 
ducing the variable y = n/n eq we obtain 
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Zcdt 



ydy 



y - 



V 



(5.5) 



and the characteristic time for this collapse process is 
thus 



L 

T PS — T£ . 



Co 



(5.6) 



A similar scenario occurs in the bridge-stretching 
regime (Fig. [l]). In this stage the p randomly placed 
pearls are roughly stationary, and they grow by assim- 
ilating monomers from the bridges joining them. The 
bridges are initially undeformed. This stage ends when 
the bridges consist of a stretched string of collapse blobs. 
The average length of each stretched bridge is R p w 
(n/g c p)£ c - As, within this model, the p bridges are ran- 
domly placed, the overall span of the chain is R ~ p x l 2 R v . 
By assuming that at this stage the overall size of the chain 
remains R ps Rf ~ N 3 ^ 5 a and that n ps N we are led to 
the result 



p : 



iV 4/5 /5c 



(5.7) 



This amounts to the completion of all short length-scale 
process prior to the onset of the overall shrinking of the 
chain. The growth of the dissipation with the length- 
scale involved justifies this estimate. Clearly, this is a 
rough estimate and, in reality, R is expec ted to shrink 
somewhat. The estimated value of p, Eq. (5.7), specifies 
R p at the end of the stretching stage to be 

R P « A^ 1/5 £ c (5.8) 

and the corresponding characteristic time is than 

r bs ps (5.9) 

In contrast with t p , the time scale t& s , depends on both 
p and N. 

VI. COLLAPSE OF THE "PEARL NECKLACE" 

At the end of the bridge-stretching regime, the config- 
uration of the chain is that of a random walk of p steps 
of length R p and with an overall radius Rf ~ N 3 ^ 5 . The 
p ps TV 4 / 5 / g c elementary steps are dumb-bells comprising 
of two pearls joined by a stretched string of £ c -blobs. It is 
impossible to incorporate more monomers into the pearls 
while keeping them stationary. The pearls can only grow 
by shortening the bridges, thus causing the overall shrink- 
ing of the chain from R ~ Rp to R ps R c ~ A 1 / 3 . It is 
assumed that the pearls are the vertices of a random walk 
and that their number, p, is constant. This approxima- 
tion is reasonable in the early phase of this stage but 
may be questionable at later times. This issue will be 
discussed later. On short time-scales the bridges equili- 
brate rapidly whereas the pearls remain essentially fixed 



in space. Consequently, the positioning of the pearls does 
not give rise to an entropy and the free energy corre- 
sponding to this process is F ps pF p b, where F p b is the 
free energy of a pearl attached to a stretched string of 
collapse blobs. In turn, F p b comprises two terms. One is 
the excess free energy of a "free" collapse blob as given 
by the kT per blob Ansatz. The second is the elastic 
free energy of the stretched bridge of length Rp with re- 
spect to its reference state, i.e., a random walk of col- 
lapse blobs. The span of an undeformed bridge with an 
ideal coil configuration is R p q ps (n/pgc) 1 ^ 2 ^ and the 
associated Gaussian stretching penalty is kT(R p /R p o) 2 . 
Assuming that the overall configuration is that of a ran- 
dom walk, R ps p 1 / 2 i? p , the free energy of the chain as a 
whole is given by 



F 



chain 

kT 



P 



N — 
P9c 



R 2 



(6.1) 



The equilibrium state of the chain, at a given instant, 
is specified by dF/dn = for constant R and p, thus 
leading to 

(6.2) 



We assumed that the decrease in R is slow in comparison 
to the relaxation of n for a given R. Accordingly, the free 
ener gy for a certain R, as obtained from Eqs. ( |6.l[ ) and 
(|6.2D is given by 



F, 



chain 

kT 



N R 

1- Vpt> 

9c Kc 



and the corresponding force, / /kT ps dF/dR, is 
kT _ 



(6.3) 



(6.4) 



where f co ps kT/^ c is the equilibrium tension in the 
dumb-bell shaped, coexistence configuration. If the 
string of pearls is perfectly collinear, R ss pR p , similar 
considerations lead to R ss (nj g c )£,c and / ps f co . In this 
situation only the terminal pearls experience a net force. 
As the tension in the bridges, ignoring finite size effects, 
is f co , the net force on the middle pearls is zero. How- 
ever, if the pearls are the vertices of a random walk, each 
pear l is subject to a random force ~ f co thus leading to 
( |6.4[ ). The thermodynamic driving force for the collapse 
is simply dF c h a in/dt ps fdR/dt. This reflects the rate of 
change in the number of "free" £ c -blobs, d(n/g c )/dt. As 
before, two dissipation mechanisms are involved. The 
hydrodynamic dissipation is T(dSh/dt) ps r](dR/ dt) 2 R. 
The dissipation modes in pearl necklace situation differs 
from the stretched-bridges case in two respects. First, 
R(t) replaces L w const' because the coil is assumed 
to shrink as a sphere. Second, as the pearls are no 
longer stationary it is necessary to consider two dissi- 
pation modes associated with the Stokes drag force. One 
is due, as before, to the contribution of the stretched 
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bridges whereas the second arises from the motion of 
the hydrodynamically impenetrable pearls. The pearls 
term is T{dS p /dt) w prjr p (dR/dt) 2 , where the volume 
of a pearl of closed packed £ c -blobs, r 3 ss (N — n)/pg c , 
specifies the radius of an individual pearl. The bridges 
term is, similarly, T(dSb/dt) rj prjR p (dR/dt) 2 1 where 
R p rj {n/g c )^ c is the length of an individual bridge. Al- 
together, dF cham /dt = -T(dS t otai/dt) leads to 



dt 
P9c — 



VP9c 



N-n 
P9c 



1/3 



n 

9c 



dn, (6.5) 



where the three terms on the right hand side correspond, 
respectively, to the dissipation due to hydrodynamic flow, 
the pearls and the bridges. In this case it is simpler to 
solve individually each dF c h a in/dt = —T(dSi/dt), in or- 
der to find the characteristic time for the different dis- 
sipation modes and then to identify the dominant one. 
The hydrodynamic dissipation yields 



Th 



P~ 3/ \, 



whereas the bridges contribution leads to 



N > _! 
n ~ 1 — ) p T£, 



and the pearls dissipation to 



Tpd 



9c J 



(6.6) 



(6.7) 



(6.8) 



Estimating pbyps N 4 ^ 5 /g c yields r h ~ N^ 5 /gl /2 , 
r b ~ N 6 / 5 /g c , and T pd ~ N 16 / 15 /g c leading to the iden- 
tification of T;, rj (iV 6 / 5 / g c )r^ as the longest relaxation 
time and the bridges dissipation as the rate-determining 
dissipation mode. Note, however, that this is the case 
only for sufficiently deep quenches, when g c < N 2 / 5 . 
By contrast, hydrodynamic dissipation is dominant when 
g c > TV 2 / 5 . 

The discussion presented above is based on three as- 
sumptions: (i) Finite-size corrections are negligible and 
f co w kT/£ c . (ii) The pearls are placed at the vertices of 
a random walk (iii) p > 1 and p ~ const' . The first as- 
sumption is easiest to justify. The finite size correction to 
f co is l/r p , where r p is the radius of the pearl. This cor- 
rection is expected to be small at the end of the bridge- 
stretching stage. In any case, it decreases with time. 
The second assumption is reasonable, within our model, 
at the onset of the collapse of the pearl necklace. Its 
range of validity is however uncertain. Some simulation 
results support this picture whilst others sug- 
gest that the chain may approach a linear configuration. 
The assumption that p ^> I is reasonable for N ^> I. 
Both [^5[ and jlTj support this view. The assumption 
of p ri const 1 is clearly a rough approximation and p is 



expected to decrease with time. It is nevertheless inter- 
esting to note that both [|l5| and JlTj suggest that p varies 
very weakly over a wide range of times. "Pearl-pearl" co- 
alescence is the most likely mechanism for the decrease in 
p. Initially, this is due to one-dimensional diffusion of the 
pearls that coalesce upon contact because of the surface 
free energy. At this point it is useful to consider two pos- 
sible scenarios for the coalescence of the pearls. In the 
first, pearl diffusion is fast in comparison to the overall 
collapse of the chain. In this case the one dimensional 
diffusion incurs no frcc-cncrgy penalties as long as n is 
constant. The favored mode thus involves the transfer 
of monomers between two adjacent bridges and through 
the adjoining pearl. This transfer of monomers across the 
pearl must take place along the trajectory of the chain. 
It is, thus, a reptation type motion in the sense that the 
monomer needs to diffuse across a distance (r p /a) 3 a in- 
stead of r p . This dissipation mode may well dominate 
over the Stokes drag force. Only the terminal pearls may 
diffuse in a non-reptative mode. In the other limit the 
pearl diffusion is relatively fast. In this case it can result 
from asymmetry in the rate of assimilating monomers 
from the two bridges. In this situation, the diffusion co- 
efficient will reflect the Stokes drag force on the pearls. 
In any case, the dynamics of diffusion-controlled coales- 
cence in one-dimensional is complex. The case of the 
A + A — ► A reaction, with constant diffusion coefficient 
D, was studied extensively [34], and a rigorous solution 
was reported only recently [35|. At short times the lin- 
ear concentration, c, follows the mean-field rate equation 
dc/dt oc — c 2 . However, the asymptotic solution at large 



times is c ~ 1/yDt instead of c ~ 1/Dt, as predicted 
by the mean-field treatment. The situation encountered 
in the present case is more complicated in two respects: 
First, the diffusion coefficient is no longer a constant. 
Rather, it depends on r p and thus on t. Second, the to- 
tal length varies with p and r p thus introducing a novel 
mechanism for the t dependence of c. A complete anal- 
ysis of this problem is beyond the scope of the present 
paper. It is however possible to estimate a characteris- 
tic time, t pp , for the onset of coalescence. t pp may be 
identified with the time required by a pearl in order to 
diffuse across a distance R p « -/V 1//5 £ c , as specified by the 
Einstein diffusion relation D p t pp ri R 2 . A lower bound 
for t pp is obtained by assuming that the diffusion coef- 



ficient of the pearl scales as D p 



l/£ c thus leading to 



VII. CONCLUDING REMARKS 

The literature concerning the theory of the dynamics of 
collapse describes three different approaches: Langevin- 
type theories, simulation studies and phenomenological 
models. It is simplest to compare the present model to 
the earlier phenomenological models. The kinetics of col- 
lapse were first considered by de Gennes in 1985 |0. In 
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this model the collapsing chain was viewed as a "sausage" 
of collapse blobs. With time, the sausage diameter in- 
creases and its length decreases. The driving free en- 
ergy was identified with the interfacial free energy and 
the dissipation was attributed to hydrodynamic flow. A 
different view was later proposed by Buguin, Brochard- 
Wyart and de Gennes jl2). Within this model the driv- 
ing free energy is the transfer free energy into the dense 
phase. Two scenarios were considered for the early col- 
lapse. In one, the pearls form the vertices of a random 
walk. This scenario differs from the present model in 
two respects. First, the bridges joining the pearls are 
assumed to be unstretched. Second, the associated dis- 
sipation arises only from the hydrodynamic flow. Im- 
portantly, the collapse time within this picture does not 
depend on the number of pearls, p. In the second sce- 
nario, the chain is envisioned as a dumb-bell consisting 
of two pearls, p = 2, joined by a stretched string of col- 
lapse blobs. The dissipation in this case is attributed to 
the Stokes drag force on the pearls. Most recently the 
problem was re-examined by Klushin Jig} ]. In this model 
the driving free energy is the transfer free energy to the 
dense phase and the dissipation is due to the Stokes drag 
force. The chain is envisioned as a self-avoiding random 
walk of Kuhn lengths. In turn, these can comprise ei- 
ther stretched bridges joining pearls or by a sequence of 
collinear pearls connected by stretched bridges. Thus, 
the number of Kuhn lengths can change either because 
of coalescence or because of the formation of a collinear 
sequence. The rate of change in the number of Kuhn 
lengths is assumed to follow unimolecular kinetics. The 
characteristic time for this process is identified with that 
of a shrinking dumb-bell, assuming a constant size of the 
pearls. The number of pearls within a collinear sequence 
is assumed to decrease because monomer transfer or coa- 
lescence, both driven by the surface free-energy. Within 
this model, the collapse is largely due to the decrease in 
p. The models of Buguin et al. and of Klushin focus, in 
effect, on the kinetics of the collapse of a pearl necklace 
and do not consider possible earlier events. 

In principle, the Langevin equation can provide a com- 
plete description of the collapse on all time-scales and 
thus allows for a systematic analysis of the problem. Un- 
fortunately, a rigorous and complete solution of the ap- 
propriate Langevin equation is difficult. Accordingly, all 
the Langevin-type theories introduce certain approxima- 
tions so as to reach a mathematically tractable formal- 
ism. The treatment is mathematically demanding, and 
often requires numerical calculations. A detailed compar- 
ison with phenomenological theories is difficult because 
of the very different formulations. It is interesting how- 
ever to note that both the theory of Garel at al. |1C]] and 
of Kuzentsov et al. Q distinguish between two regimes 
in the early collapse: first, a very fast process, whose 
characteristic time is independent of N and is thus in- 
terpreted as local rearrangement; and second, a much 
slower process, reflecting large scale reorganization, with 
a characteristic time that increases with N. 



We hope that a full picture of the kinetics of col- 
lapse will eventually emerge from molecular dynamics 
simulations involving a polymer chain as well as solvent 
molecules. The currently available simulation studies do 
not account fully for hydrodynamics thus introducing an 
uncertainty in the interpretation of the results. In cer- 
tain studies the interpretation is further hampered by the 
imposition of an "all trans" initial state p^ , |l6| . With 
this caveat, the simulations mostly support the notion 
that the collapse involves a string of pearls joined by 
stretched bridges as an intermediate stage. The configu- 
ration adopted by the string of pearls is less clear. Snap- 
shots of the chain configuration as obtained by Langevin- 
equation simulation of Byrne et al. Q suggest an almost 
collinear string of pearls. The visualization of configura- 
tions obtained by Monte Carlo simulation of Kuznetsov 
et al. Jljj reveal a random walk pearl necklace. These 
last two studies also suggest that the first step of the 
collapse involves the formation of localized clusters, and 
that subsequent cluster growth is attained by assimilat- 
ing monomers from the bridges thus causing stretching. 
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